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ABSTRACT 

The HH 24 MMS protostellar system was observed in the 6.9 mm continuum with a high angular 
resolution (0.5"). HH 24 MMS was resolved into two sources. The separation between sources 1 and 2 
is —0.9" or 360 AU. The spectral energy distribution suggests that the 6.9 mm flux is almost entirely 
from dust. The 6.9 mm image and the spectrum suggest that HH 24 MMS may be a protostellar binary 
system. Total mass including the accretion disks and the inner protostellar envelope is ~1.4 Mq. Disk 
masses of sources 1 and 2 are 0.6 and 0.3 Mq, respectively. Both sources are highly elongated. The 
difference in the position angle of the two disks is —45°, which means that HH 24 MMS is a highly 
misaligned protobinary system. The misalignment suggests that turbulent fragmentation may be the 
formation mechanism relevant to the binary systems with a separation of a few hundreds of AU, such 
as the HH 24 MMS system. 

Subject headings: accretion disks — binaries: general — ISM: individual (HH 24 MMS) — ISM: 
structure — stars: formation 



1. INTRODUCTION 

Young stellar objects are often classified using the spec- 
tral energy distribution (SED; Lada 1987; Andre et al. 
1993). Class sources have cold spectra that peak in the 
submillimcter/far-IR band. Class sources are inter- 
preted as highly embedded protostars still in their main 
accretion phase (Andre & Montmerle 1994). A protostar 
is thought to contain a star-disk system surrounded by a 
collapsing envelope (e.g., Shu et al. 1987). It is impor- 
tant to understand the structure of Class objects since 
they represent the earliest phase of star formation. How- 
ever, it is not easy to study such systems owing to the 
high extinction in the optical or IR wavelength shorter 
than ~20 fim. Observations in the millimeter band pro- 
vide us images of protostellar envelopes and deeply em- 
bedded accretion disks. High-resolution millimeter im- 
ages often show that many Class sources are multiple 
protostellar systems (Looney et al. 2000). 

Many stars form as members of binaries or multiple 
systems (e.g., Duquennoy & Mayor 1991). The multi- 
plicity of young pre-main-sequence stars is higher than or 
comparable to that of main-sequence stars (Duchene et 
al. 2007), implying that multiple systems form through 
fragmentation of cores in the early stage of star forma- 
tion. Several kinds of fragmentation mechanisms have 
been suggested, including rotational fragmentation, tur- 
bulent fragmentation, and disk fragmentation (Goodwin 
et al. 2007). Fragmentation of cloud cores has been in- 
vestigated through numerical simulations in various sit- 
uations including many kinds of physical and chemical 
effects, and it was found that most cores tend to frag- 
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ment into multiple objects (see Goodwin et al. 2007 and 
references therein) . While collapsing cores fragment eas- 
ily, the relative importance of rotation and turbulence in 
star forming regions remains unresolved. Observations of 
Class sources in the millimeter band can provide clues 
to the fragmentation mechanism. 

HH 24 MMS is a Class source located -20" south- 
west of the Herbig-Haro object HH 24A in the L1630 
cloud (Chini et al. 1993; Ward-Thompson et al. 1995; 
Bontemps et al. 1996). It was detected in the 3.6 cm and 
2.0 cm continuum by Bontemps et al. (1996). Chandler 
et al. (1995) suggested that the dust emission of HH 
24 MMS originates from two components: an unresolved 
disk and an extended envelope. Bontemps et al. (1995) 
suggested that the 3.6 cm emission at the position of 
the submillimeter peak originated from thermal free-free 
emission, but later they found that the unresolved com- 
ponent of the centimeter continuum source has an SED 
consistent with the dust disk (Bontemps et al. 1996). 
Reipurth et al. (2002) imaged the 3.6 cm source with an 
angular resolution of —1". High-density molecular gas 
was detected in molecular lines (Gibb & Heaton 1993). 
Non-detection of X-rays from HH 24 MMS suggests that 
either HH 24 MMS may be very deeply embedded or HH 
24 MMS does not have (proto)stellar flares to produce 
X-ray emission (Ozawa et al. 1999; Simon et al. 2004). 

Despite the name, most Herbig-Haro objects in the HH 
24 complex are unrelated to HH 24 MMS but belong to 
outflows driven by SSV 63 and its companions (Mundt 
et al. 1991; Davis et al. 1997). Near-IR observations 
revealed that HH 24 MMS is driving its own outflow or 
an "H2 jet" oriented in the northeast-southwest direction 
(Bontemps et al. 1996; Davis et al. 1997). Eisloffel 
& Mundt (1997) suggested that HH 24A and HH 24L 
may belong to the HH 24 MMS jet. This jet was also 
detected in the radio continuum and high-velocity CO 
gas (Bontemps et al. 1996). 

In this paper, we present high resolution images of the 
Class protostar HH 24 MMS observed in the 6.9 mm 
continuum using the Very Large Array (VLA) . In § 2 we 
describe our observations. In § 3 we report the results 



2 



Kang et al. 



TABLE 1 

Continuum Sources in the HH 24 MMS Region 





Peak Position 


Peak Intensity 


Source SizE a 




Total FLux a 


Mass 


Source 


Q!2000 


#2000 


(mjy beam -1 ) 


FWHM 


P.A. 


(mJy) 




1 

2 

3 


05 46 08.37 
05 46 08.42 
05 46 08.56 


-00 10 43.7 
-00 10 43.1 
-00 10 42.9 


2.54 ± 0.07 
0.87 ± 0.07 
0.45 ± 0.07 


0.72 ± 0.04 x 0.28 ± 0.03 
1.35 ± 0.23 x 0.18 ± 0.05 
unresolved 


-66 ± 3 
69 ± 2 


4.9 ± 0.4 
3.0 ± 0.6 


6+ - 6 
u - D -0.3 

3 +0 - 4 



Note. — Units of right ascension arc hours, minutes, and seconds, and units of declination arc degrees, arcminutes, and arcseconds. Source 
parameters are from the C-array uniform-weight map (Fig. lb). Fluxes were corrected for the primary beam response. 

a Parameters of the best-fit two-component elliptical Gaussian profile. Deconvolved size is in arcseconds, and position angle is in degrees. At 400 
pc, the deconvolved sizes of sources 1 and 2 correspond to 290 AU X 110 AU and 540 AU X 70 AU, respectively. 



of the 6.9 mm imaging. In § 4 we discuss the physical 
properties of the dust continuum sources in the HH 24 
MMS region. 

2. OBSERVATIONS 

The HH 24 MMS region was observed using the VLA of 
the National Radio Astronomy Observatory 6 in the stan- 
dard Q-band continuum mode (43.3 GHz or A = 6.9 mm) 
in two observing tracks. The phase and amplitude were 
determined by observing the nearby quasar 0532+075 
(PMN J0532+0732). The phase tracking center was 
c*20oo = 05 h 46 m 08.38 s and 6 20 oo = -00°10'43.36". The 
first track was carried out in the D-array configuration 
with twenty- five antennas on 2003 March 30. The flux 
density of 0532+075 was bootstrapped from the quasar 
0542+498 (3C 147) assuming that its flux density is 0.72 
Jy, which is the flux density measured within a day of our 
observations (VLA Calibrator Flux Density Database 7 ). 
The bootstrapped flux density of 0532+075 is 1.05 Jy. 
The second track was in the C-array configuration with 
twenty-three antennas on 2004 March 2. The flux cal- 
ibration was done by observing the quasar 0713+438 
(QSO B0710+439). The flux density of 0713+438 was 
set to 0.20 Jy, which is the flux measured within 2 days 
of our observations according to the VLA Calibrator Flux 
Density Database, and the bootstrapped flux density of 
0532+075 is 0.70 Jy. Maps were made using a CLEAN 
algorithm. 

3. RESULTS 

Figure la shows the map of the 6.9 mm continuum 
toward the HH 24 MMS region from the C- and the D- 
array data combined. The peak intensity is 5.66 mJy 
beam" 1 . HH 24 MMS was clearly detected and shows 
an extended structure. 

Detailed structures of compact objects were revealed 
by imaging with the C-array data only (Fig. lb). HH 24 
MMS was resolved into two objects, and another source 
was detected ~3" east of HH 24 MMS. Their coordinates 
and peak intensities are listed in Table 1. Source 1, the 
brightest one, is near the 3.6 cm continuum source de- 
tected by Reipurth et al. (2002) and is elongated in the 
southeast- northwest direction. In contrast, source 2 is 
elongated in the northeast-southwest direction. The sep- 
aration between sources 1 and 2 is ~0.9" or 360 AU at a 
distance of 400 pc (Anthony-Twarog 1982). 

6 The NRAO is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, 
Inc. 

7 See http://aips2.nrao.edu/vla/calflux.html 



Source 3 was detected over the detection limit of 0.4 
mJy beam" 1 (S/N = 6) in the C-array map (Fig. lb), 
and it can also be seen in the C and D-array map (Fig. 
la) as a weak peak. Source 3 is unresolved. Since it was 
not detected in other wavelengths, the nature of source 
3 is unclear. 

3.1. Spectral Energy Distribution and Mass 

In low-mass star forming regions, continuum emission 
in the centimeter-millimeter band is usually either free- 
free radiation from hot ionized gas or thermal radiation 
from cold dust, or a combination of both. The free- 
free emission usually has a small (<1) spectral index, 
while the dust emission has a large (>2) spectral index 
(Reynolds 1986; Anglada et al. 1998). The SED was ex- 
amined to investigate the emission mechanism of HH 24 
MMS. Since sources 1 and 2 are resolved in the 6.9 mm 
band only, SED of the whole HH 24 MMS system was 
constructed. The total flux of the HH 24 MMS system, 
measured in a 4" x 4" box, is 12.0 ± 0.5 mJy (from the 
map shown in Fig. la), and the SED in the wavelength 
range of 3.6 cm to 3.4 mm is shown in Figure 2. 

We first tried to fit the SED of HH 24 MMS as a sum of 
two power-law components, one from the free-free emis- 
sion and the other from the dust emission, over the wave- 
length range of 3.4 mm to 3.6 cm (8.3-88 GHz). The best 
fit has spectral indices of 2.1 ± 0.3 and 3.7 ± 0.2 (Fig. 
2 left panel). It was found that acceptable fits can only 
have large (>2) spectral indices for both components, 
which clearly means that the emission is mostly from 
dust and that the free-free emission is negligible in the 
wavelength range considered. 

Then, as the dust emission seems to be dominating 
even in the centimeter wavelengths, the SED was fitted 
again with a single-component power-law distribution. 
The fits were made using three data points, excluding 
the 3.6 cm data, to avoid any contribution from free- free 
emission. The best fit has a spectral index of a = 2.8 ± 
0.2 in the millimeter range (Fig. 2 right panel). 

Extrapolating the best-fit SED to longer wavelengths, 
a substantial fraction (~90%) of the 3.6 cm flux seems 
to come from the dust emission. Therefore, the free-free 
emission, if any, must be very weak. If 10% of the 3.6 
cm flux comes from optically thin free-free emission, its 
contribution to the 6.9 mm flux would be only ~0.2% (or 
0.02 mJy). That is, the 6.9 mm flux of HH 24 MMS is 
almost entirely from dust. 

As the millimeter continuum flux is mostly from dust, 
the mass of the emission structure can be estimated from 
the SED. To derive the mass from the dust continuum 
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Fig. 1. — Maps of the A = 6.9 mm continuum toward the HH 24 MMS region, (a) Natural-weight CLEAN map from the C- and the 
D-array data combined. The contour levels arc 1, 2, 4, 8, and 16 times 0.24 mjy beam -1 . The rms noise is 0.08 mjy beam -1 . Dashed 
contours represent negative levels. Shown at the bottom right corner is the synthesized beam: FWHM = 1.0" X 0.8" and P.A. = -6°. 
The straight line at the bottom left corner corresponds to 500 AU at a distance of 400 pc. (6) Uniform-weight CLEAN map from the 
C-array data only. The contour levels are 1, 2, 4, and 8 times 0.21 mjy beam -1 . The rms noise is 0.07 mjy beam -1 . Source numbers are 
labeled. The synthesized beam has FWHM = 0.5" X 0.4" and P.A. = -39°. Plus sign: Position of the 3.6 cm continuum source (Rcipurth 
et al. 2002). The size of the marker corresponds to the synthesized beam, (c) Best-fit image of the complex containing sources 1 and 2 
from the two-component elliptical Gaussian fit to the map shown in panel (6). Thick ellipses: Deconvolved size (FWHM) of each Gaussian 
component. 



flux, the mass opacity given by Beckwith & Sargent 
(1991) is assumed, 



k v = 0.1 ( — ] cm z 



(1) 



where v is the frequency, vq = 1200 GHz, and (3 is 
the opacity index. This opacity is for total mass of 
gas and dust. The opacity index of HH 24 MMS is 
[3 fa a — 2 = 0. 8. This small value of [3 implies the pres- 
ence of large dust grains probably caused by the grain 
growth in the high density environment (e.g., Miyake & 
Nakagawa 1993), which was pointed out by several au- 
thors previously (Chandler et al. 1995; Ward-Thompson 
et al. 1995). 

The dust temperature, T d , used in previous works 
ranges from 12 K to 100 K (Phillips et al. 2001; Chan- 
dler et al. 1995). Here we adopt T d = 20 ± 5 K from 
Ward-Thompson et al. (1995). In the 6.9 mm maps, the 
peak brightness temperature is about 8 K, and the mean 
brightness temperature in a 4" x 4" box is less than 1 
K. Since the brightness temperature is much lower than 
the dust temperature, the optical depth of the 6.9 mm 
continuum emission must be small. Assuming optically 
thin emission, the mass can be estimated by 



M 



F V D 2 



n v B„(T d y 



(2) 



where F v is the flux density, D is the distance to the 
source, and B v is the Planck function. Assuming D = 
400 pc, the total mass is 1.4+J'y M@, which is consistent 
with the mass estimate of Chandler et al. (1995). This 
mass includes both the accretion disk(s) and the inner 
protostellar envelope. 

4. DISCUSSION 

4.1. Decomposition of the HH 24 MMS System 

At the first glance of the 6.9 mm maps, the HH 24 
MMS system appears like a combination of a protostar 



(source 1) and an outflow (source 2). However, the SED 
shows that, for both sources 1 and 2, the 6.9 mm flux is 
mostly from dust. Especially, the 6.9 mm flux of source 
2 must be entirely from dust because the 3.6 cm source 
is not associated with source 2 (Fig. lb). That is, it is 
very unlikely that source 2 is a thermal jet. 

Since both sources are extended and partially overlap- 
ping, a simple two-component fitting to the image was 
done to decompose the system. The image was fitted 
with a sum of two elliptical Gaussian distributions. Only 
the pixels around source 1 and 2 with intensities above 
0.3 mjy beam -1 were considered for the fit, and source 3 
was ignored. Table 1 lists the parameters of the best-fit 
Gaussian profiles, and Figure lc shows the best-fit inten- 
sity distribution, which agrees with the observed image 
within three times the rms noise. The total flux of the 
two-component Gaussian profile also agrees with that of 
the real image within 0.1 mjy. 

There is a bipolar jet seen in H2 emission in the HH 
24 MMS region with the axis in the northeast-southwest 
direction (Bontemps et al. 1996). The near-IR H2 map 
of Davis et al. (1997) shows that HH 24 MMS is almost 
certainly the driving source of this outflow. The position 
angle of these H 2 knots (NE 1, NE 2, and SW 2) with 
respect to HH 24 MMS is 32° ± 6°, which is perpendic- 
ular to the elongation of source 1. Therefore, source 1 
seems to be the driving source of the bipolar outflow. 

The peak position of the 3.6 cm source (Reipurth et 
al. 2002) does not exactly coincide with the position of 
source 1 (Fig. lb). The position displacement is about 
0.7". Since the 3.6 cm source is weak and extended 
(~1")> it is not clear that the displacement is significant. 
In this paper we assume that the 3.6 cm source is identi- 
cal to or closely associated with source 1, but we cannot 
rule out the possibility that the 3.6 cm source indicates 
the existence of yet another young stellar object. This 
issue should be resolved by more sensitive observations 
of the centimeter continuum in the future. 

For source 2, three kinds of interpretation are possible. 



4 



Kang et al. 




1.2 1.4 l. 
Log v [GHz] 



1.2 1.4 1.6 
Log v [GHz] 



Fig. 2. — Spectral energy distribution of HH 24 MMS (sources 1 and 2 combined). Flux densities are from Reipurth et al. (2002), 
Chandler et al. (1995), and this work. Solid curves: Best-fit power-law spectra. Left panel: Two-component power-law spectrum. The 
best-fit spectral indices are a\ = 2.1 ± 0.3 and a.2 = 3.7 ± 0.2. All the four data points were considered for the fit. Dotted lines: Each 
component of the fit. Right panel: Single-component power-law spectrum. The best-fit spectral index is a = 2.8 ± 0.2. The 8.4 GHz data 
point was not considered for the fit. 



(1) Source 2 is another protostar with a projected sepa- 
ration of 360 AU from source 1. (2) Source 2 is a clump 
of dust swept and heated by the outflow. (3) Source 
2 has multiple clumps in the dense envelope as seen in 
other protostars such as NGC 1333 IRAS 4B (Looney et 
al. 2000). The first possibility seems to be the simplest 
explanation. The second possibility is rather unlikely 
because the position angle of the elongation of source 2 
(~69°) is different from the direction of the outflow by 
~40°. In addition, the mass of source 2 seems too large 
for a clump around outflow. The third possibility can 
only be tested by imaging with a higher resolution, but it 
is rather peculiar that only one side of the envelope has 
dense clumps. Therefore, we prefer the interpretation 
that HH 24 MMS may be a protostellar binary system, 
and we will adopt this interpretation in the following dis- 
cussions. 

The separation between sources 1 and 2 is within the 
range of the expected separation of young binaries. In 
a few examples of well-studied Class protostars, obser- 
vations and modeling suggest that the centrifugal radius 
ranges from 34 AU to 600 AU (see Myers et al. 2000 and 
references therein). Reipurth et al. (2004) carried out 
a radio-continuum survey of embedded objects with dis- 
tances between 140 and 800 pc and found many binaries 
with separations ranging from 0.5" to 12". Many numer- 
ical simulations of core fragmentation produce multiple 
systems with separations of a few hundreds of AU (Good- 
win et al. 2007). 

Both components of HH 24 MMS are highly elongated, 
with axial ratios of 3 for source 1 and 8 for source 2, which 
suggests that each of them may be a disk-like flattened 
structure viewed nearly edge-on. The size of each compo- 
nent is somewhat larger than the expected size of an ac- 
cretion disk supported by Keplerian rotation (^100 AU), 
and each of the 6.9 mm continuum sources may consist 
of a compact accretion disk surrounded by a flattened 
inner envelope or a pseudodisk (a contracting disk-like 
structure that is not completely supported by rotation; 
Galli & Shu 1993; Tomisaka 2002; Allen et al. 2003)! 
Assuming that the dust properties of each source are not 



very different from the overall dust properties of the HH 
24 MMS system, the mass can be derived from the flux 
of each source (Table 1). 

4.2. Misaligned Disks 

The majority of stars belong to multiple-star sys- 
tems (Duquennoy & Mayor 1991), and understanding of 
star formation cannot be complete without knowing the 
physical processes in multiple-protostar systems (Tohline 
2002). Study of the protobinary in the earliest evolution- 
ary stage is essential in understanding how the fragmen- 
tation occurs and what is the initial condition of the bi- 
nary evolution. There are two main mechanisms for frag- 
mentation: (1) rotational fragmentation and (2) turbu- 
lent fragmentation (Goodwin et al. 2007). The outcome 
of these two mechanisms can be quite different. The for- 
mer occurs when the kinematics of the initial core is dom- 
inated by systematic rotation, and the resulting binary 
would be well-ordered, i.e., disks would be well-aligned. 
In contrast, the latter occurs when the core is dominated 
by random motions, and the resulting binary is likely to 
have misaligned disks. Therefore, the alignment of disks 
is one of the key observational measures that can provide 
strong constraints on the fragmentation mechanism. 

In principle, multiple outflows can give some implica- 
tions on the degree of binary disk alignment because out- 
flows are ejected perpendicular to the circumstellar disks. 
There are some examples of misaligned systems inferred 
from jets emanating with different position angles (Gre- 
del & Reipurth 1993; Reipurth et al. 1993). However, 
this method is problematic because two well-aligned out- 
flows would be misidentified as a single outflow. The best 
way to test the disk alignment may be the direct imaging 
of both disks. 

If HH 24 MMS is a binary system as discussed in the 
previous section, it is one of the rare examples of proto- 
binary system imaged with an angular resolution high 
enough to investigate the degree of alignment. Each 
source may consist of a compact accretion disk sur- 
rounded by a pseudodisk, and the axis of the accretion 
disk is very likely to be parallel to the axis of the pseu- 
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dodisk around it. Then the axis direction of the accretion 
disk can be inferred from the minor axis of the elliptical 
fit to the 6.9 mm source. 

The major axes of the two components of the HH 24 
MMS system are not aligned at all. The difference in 
the position angle between the binary components, from 
the Gaussian fit (Table 1), is ^45° (or ~135° if they are 
retrograde). HH 24 MMS makes a stark contrast to the 
other known example of resolved binary disks, LI 551 IRS 
5. In the L1551 IRS 5 system, the two disks are very well 
aligned (Rodriguez et al. 1998). The difference between 
the two systems is probably owing to the binary separa- 
tion: 360 AU in HH 24 MMS and 45 AU in L1551 IRS 
5. In the classification scheme of multiple young stellar 
objects suggested by Looney et al. (2000), HH 24 MMS 
belongs to the common-envelope type while L1551 IRS 
5 belongs to the common-disk type. We speculate that 
the rotational fragmentation is important in close bina- 
ries and the turbulent fragmentation is important in wide 
binaries, though more examples are certainly needed to 
draw a statistically significant conclusion. Other exam- 
ples may include IRAS 16293-2422 and NGC 1333 IRAS 
4A (Looney et al. 2000), but their structures are corn- 
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plicated and their disk orientations are less obvious. 

Though studies through direct imaging is difficult, in- 
direct evidences show that misaligned systems are not 
necessarily unusual. High-resolution near-IR polarime- 
try revealed plenty of examples of misaligned binaries in 
the Class II phase (Monin et al. 2007). Disks in T Tauri 
binaries are not perfectly coplanar (Jensen et al. 2004), 
and the relative alignment of the disks may decrease after 
formation (Bate et al. 2000). Therefore, HH 24 MMS can 
be considered as a precursor of such misaligned T Tauri 
systems. As HH 24 MMS is a rare example of a wide 
binary that can be imaged well enough to study the disk 
alignment, this system gives us a chance to examine the 
initial conditions of binary evolution. 
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